NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



^^Anmmaij (Mm^ 

American Institute of Aeronautics and Astronautics 
1290 Avenue of the Americas, New York, N. Y. 10104 
Telephone 212/581-4300 

June 29, 1981 
TO: Shirley Peigara 

FROM; Norma Brennan, Director, Editorial Department 

A back-up paper is enclosed for the following Synoptic: 

Author(s): W. H. Hui and Murray Tobak 

Title of Synoptic: “Unsteady Newton-Busemann Flow Theory Part II: Bodies of 

Revolution" (Log No. J12130) 

Title of Back-up Paper: Same as above. 



Correspondence with: Dr. Murray Tobak 

Aerothermodynamics Section 
Entry Technology Branch 
Attn: STP: 234-1 
NASA Ames Research Center 
Moffett Field, California 94035 



(NASA>TM-804i9) ONSIEADi NEhlOM-BUSE BANN M81-28J71 

fLOi THEOBY. i*ABT 2i cO£i£S Of aEVOLUTlOM 
(NASA) J7 y tiC A03/MF AUl CSCX, 20D 



UNSTEADY NEWTCW-BUSEMMIN FLOW THEORY 


PART II: BODIES OP REVOLUTION 


By 


W. H. Hui and Murray Tobak 


/April 1981/ 



Backup Document for AIAA Synoptic Scheduled 
for Publication in the AIAA Journal, October 1981 


NASA Ames Research Center 
Aero thermodynamics Section 
Entry Technology Branch 
Attn: STP: 234-1 

Moffett Field, CA 94035 
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PART ll: BODIES OF REVOLUTION 
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Abatract 

Newtonian flow theory for unsteady flow past oscillating bodies of revolu- 
tion at very hi^ Mach niodbers is cMipleted by adding a centrifugal force 
correction to the iapact pressures. Exact fomulas for the unsteady pressure 
and the stability derivatives are obtained In closed form which require only 
numerical quadratures of terms Involving the body shape. They are applicable 
to bodies of re^n^lutlcm that have arbitrary shapes, arbitrary thicknesses, and 
either sharp or blunt noses. The centrifugal force correction arising from 
the curved trajectories followed by the fluid particles In unsteady flow Is 
shown to be very liq>ortant. It cannot be neglected even for the case of a 
circular cone, for which tl» centrifugal force contribution In steady flow is 
zero. With the Inclusion of this correction the present theory is shown to be 
In excellent agreement with experimental results for sharp cones and for cones 
with small nose bluntness; however, the theory gives poor agreement with the 
results of experiments In air for bodies having moderate or large nose blunt- 
ness. Finally, the pitching motions of slender power- law bodies of revolution 
are shown to be always dynamically stable according to Newton-Busemann theory. 

Index categories: Nonsteady Aerodynamics; Supersonic and Hypersonic 

Flow; Entry Vehicle Dynamics and Control. 

*Vlslting Professor. Permanently Professor of Applied Mathematics, 
University of Waterloo, Waterloo, Ontario, Canada. 
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1. Introduction 



This paper continues the study of unsteady Newton-Busemsnn flou theory 
begun by the authors In Part 1 (Ref. 1) by extending it to cover the case of 
bodies of revolution as well as airfoils. Speciflcallyt we consider high Mach 
nuid>er flow past a body of revolution of arbitrary thickness and shape having 
a sharp or blunt nose, i^ich is perfomlng a c«Bbined pitching and plunging 
oscillation of low frequency. The sin is to calculate the restilting unsteady 
flow on the basis of the Newton>Bu8eaann flow aodel. In this model we conceive 
of the surface pressure as being composed of two parts: the Newtonian impact 

pressure and a centrifugal-force correction owing to the curved trajectories 
that fluid particles follow along the surface subsequent to their lnq>act. 

In Sec. II the Newtonian impact pressure is derived, and in Sec. Ill the tra- 
jectories of the fluid particles are determined. These are iised in Sec. W 
to find the steady and unsteady centrifugal force corrections. In Sec. V 
these general results are applied to study the dynamic stability of classes 
of bodies of revolution, incl\iding the effects of the unsteady centrifugal 
force and of nose bluntness. Comparisons with experimental results are also 
shown, and, finally, conclusions arc presented in Sec. VI. 

II. Newtonian Impact Pressure 

Consider a body of revolution in a uniform hypersonic flow U„ which 
performs a harmonic pitching oscillation with zero mean angle of attack about 
the pivot axis C. Simultaneously, the pivot axis undergoes a harmonic plung- 
ing oscillation, so that the body is in a coid>lned pitching and plunging 
motion (Fig. 1). The body may have a sharp or blunt nose, but cannot be 
flat-nosed. Let a fixed system of Cartesian coordinates OXTZ (with corre- 
sponding unit vectors denoted by and be such that 0 is at the 

mean position of the nose and OX is along the direction of the free stream 
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Um, Let the body-fixed eyetee of Certeelan coordinates O'xyz (with corre- 
epondlng unit vectors denoted by It end Ic) be such that C' la at the 
nose and O'x ia along the axis of syssMtry of the body. Concurrently, a 
body-fixed systM of cylindrical coordinates (x.r.4) (with corresponding unit 

, Sp, and e^) will be used which is related to (x«y.a) by 
X ■ X , y ■ r cos ^ , z » r sin f (1) 

In the following, all the lengths X, Y, Z, x, y, z, r, h, etc., are 
scaled by the length of the bndv i, velocities v by U«, density p by 
p«, pressure p by P.U^, uuU the tine variable t by i/Uo,. Thus, for 
example, P«, “ ■ £ ■ 1. 

The harmonic pitching oscillation of the body with frequency u may be 
represented by the displacement angle of pitch 6 (Fig. 1) 

0(t) - e (2) 

%ihere 


ui£ 

U.. 


(3) 


is the reduced frequency and 6 is the amplitude of oscillation. Likewise, 
the harmr'ilc plunging oscillation of the body at the same frequency may 
be represented by the linear displacOMnt (Fig. 1) where 

ikt 


Yc<t) - Tfc • 


which is related to the fllghtpath angle y(t) by 

• .. j.x ikt _ - ikt 

Yj.(t) - IkYj, e s Y(t) -ye 


(4) 


(5) 


It will be ass\sed that |6|, |y|* and k are all «1, and all terms of 
0(6^, Y^» ^Y* yk^) and higher will be neglected. 

The unit vectors of the coordinate systois Introduced above are obviously 
related by 

t - 1 - ej . J - el + t , i-t (6) 
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co« ♦ + ^ •!» ^ ■ “1 *10 ♦ + tt co« ^ 

Also, due to the notion of the body we have 


dt ^ 
It " 





0 


(7) 


( 8 ) 


Equations (6) -(8) will be used freely throughout the paper. Let the equation 
of the surface of the body of revolution be given by 

r - f(x) (9) 

where, evidently, f(0) - 0 in the chosen coordinate systos. Then the unit 
nomal vector n on the surface is 

n - y(x)[e^ - lf*(x)l (lOa) 


or 

n ■ y(x)[-lf'(x) +3^ cos ♦ + It sin ♦) (lOb) 

where 

y(x) - [I + f'*(x)r‘^* (11) 

The position vector of a fluid particle P on the body surface can be 
written as 

p - tx + JTY + (12) 

idiere 

X - X - 9f(x)cos ♦ (13a) 

Y - Y^(t) + e(x - h) + f(x)cos ♦ (13b) 

Z ■ f(x)sin ^ (13c) 

In order to obtain the tlevtonlan lapect pressure we need to find the normal 
component of the particle velocity iMtediately before and after impact. The 
nomal c<^ponent of velocity before impact Is 

(Vq)^^^^^^ - 1 • n • -y(x)lf'(x) + 0 cos (14) 
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After impact, a fluid particle movea along the surface and the normal com- 
ponent of ita velocity C^n^^fter ■ point (x,#) is equal to that of the 
body surface velocity ^ itself. Since 

- Jy + S[-lf(x)cos ♦ + J(x - h)] (15) 

fixed 



we get 

^'^n^after " ♦(5[x- h+ f (x)f ' (x)] -I- y> 

Accordingly, the Newtonian Impact pressure is 

■ fVb,£ore‘<''n’.£t.J* * + 2 f<*)c 08 tleCt) 

+ Y(t)l + 2f'(x)coi « 9(t)[x - h + f(x)f(x))) (17) 

III. Particle Trajectories 

As remarked in Ref. 1 (Part I of this study), the Ne%rtonlan impact pres- 
sure represents the pressure at the outer edge of the Newtonian shock layer, 
but the pressure at the body surface must contain an additional centrifugal- 
force correction owing to the curved trajectories followed by the fluid par- 
ticles sxibsequent to their Impact. In this section we determine tlw particle 
t.'ajectorles due to the unsteady motion of the body based on the Ne%rtonian 
assumption that the tangential cogq>onents of acceleration are zero. We choose 
two mutually orthogonal unit tangential vectors and Tj as follows: 

tj - u(x)[l + ejf'(x)l ■ u(x)[t + Jf'(x)cos ^ + lcf'(x)ein (I8a) 

Tj - e^ ■ -J sin ♦ + ^ cos ♦ (I8b) 

The velocity of a particle may be obtained by taking a tiiK derivative of 
Bq. (12) with the use of Eq. (13). We have 

V ■ Tii - df(x)co8 ♦] + J[y + 0(x - h) + x£'(x)cos ^ 

- ^f(x)sin + ^[j^f'(x)sln ^ 4- ^f(x)cos ♦) (19a) 
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or 


« Tjji -f SwtCx - ' flcoo i ^ Yvf coi ♦! + Tjl-OCx - h)ain ^ 

- Y aln ♦ + ♦£]+ n[6vi(££' + x - h)cos ♦ YW coa 4] 


(19b) 


The velocity u of a particle relative to the body auqr be obtained from 
Eq, (19b) with 6 ■ Y “ 0. Thua, 

u - Tx + Tx^f (x' (20) 

The acceleration a of a particle may then be obtained by taking a time 

derivative of velocity, yielding 

a - T.a^ + T.a + na^j (21) 

1 Tx » Tj 


where 


i - 4 - - w4^ff ' 20M4f Bin 4 + Yuf ’ cos 4 

Tj at \y/ 


*T " T ^ “ (y + 20x)ain 4 

- y|i*f" - 4*f + 26 coa 4 - 4ff sin 4)| + YW cos 4 


4^j + Yi 


(22a) 

(22b) 

(22c) 


The equatlona governing the motion of the particle are thua 
^ 264wf sin 4 + YMf * cos 4-0 


^ (f*4) - (y + 20x)f sin 4 - 0 


(23) 

(24) 


The initial condition required to determine the particle motion la Its velocity 
immediately after impact. Before impact, a particle's velocity la t: 

I - Xivil - 0f cos 4) + TjO sin 4 -ny(f’ + 0 cos 4 ) (25) 

After impact, the ps'tlcle loses the normal component of Its velocity; hence 
its velocity Vj^ Immediately after li^pact at (x^, 4i» ti) Is 

" ^iW(x^)[l - e(t^)f*(x^)cos 4^1 + T 20 (t^)sin 4 ^ (26) 
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Comparing Eq. (26) with the general expression Eq. (I9b) for velocity of a 
particle in the tangent plane, v - nv • n, me get the following initial condl 


tionst at X - x^, ♦ ■ t ■ 

k - {(e + ‘ fjDcos (27a) 

if - [(9 + y) + S(x^ - h)]sin (27b) 

In the special case when 8 - y “ 0, the flow is axisyao»tric and we get 
^ For small 6 and Yt we have 

♦ - + 0(8.y) (28) 

We therefore look for solutions to Eqs. (23) , (24) , and (27) in the form 

^ - A(x) + [9E(x) + §C(x) + yD(x) + yE(x)]cos (29a) 

♦f - [9F(x) + k(x) + yH(x) + YK(x)]sin (29b) 

Substituting Eq. (29) into Eqs. (23), (24), and (27). using the relation 

— - Q’ (x)x • uQ'{A + (0B + 0C + yD + yE]co8 (> . ) (30) 

dt 1 


and equating like terms in 0, fi, y. Y. etc. , we obtain the following set of 
ordinary differential equations: 


AA' - 0 ; 

A(xi) - m 

(31a) 

(AB)' - 0 ; 

B(x^) - -Uifi 

(31b) 

(AC) ' + ® - 0 ; 

C(xi) - -Ui[(Xi - h)fi - fil 

(31c) 

(AD)' - 0 ; 

D(x^) • -yifi 

Old) 

+ - + f ' - 0 ; 
V 

E(xi) - 0 

Ole) 


X ■ 


^Subscript 
x^» e. g. , f 


i will denote the value 

, - df(x) 
i = dx 


X«Xj^ 


of a function evaluated at 
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0 ; 

r(xp - 1 

(310 

0 } 

G(xt) ■ Xi ' h 

(31r) 

0 ; 

H(xi) - 1 

(3lh) 

0 ; 

K(x^) - 0 

(3li) 


The eolations to Eqs. (31) ere 
A • a(x£) 

B • D > -v(x^)f'(xi) 

C(x,xj^) - v>(x^){f(X|^) - (x^ - h)f'(xi)] + £* (xj[^)So(x,Xi) 


E(x.xj^) 


F(x,x^) 


G(x,xj^) 


f(x) - £(x<) 

T(5iJ5 f'(xt)s„<x.xi) 

f(Xj) 

H(x,x^) - 


ik) 


K(x,x^) - 


r f(Xi) 1 

T) [(*1 - ' Jx^ 

^ ^x £(C) - f(xp 

cjrn^ Jx^ 


where 


Jn(x,xi> -J* 


, n - 0,2 


x^ V(0£ (O 
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IV. Centrifugal Force Correction 
Aa In Ref. I, the comrifngal preaaure Is given by 

/ . It* 

a (itln ( tA I 

•n ' 

P 

vhere p la density, P is the point on the bjdy surface where pressure is 
evaluated, and P* is the point at the outer edge of the shock layer across 
froB P. While an expression for Sg, the nonal coaponent of acceleration, 
is readily found froa Eq. (22c) plus Bqa. (29) and (32), we need an expression 
for pdn. This is obtained frott the law of conservation of aasa. In a frame 
of reference fixed with respect to the body, the continuity equation reads 

^ (pdn d#) + V • (u pdn d9) - 0 (35) 

where dn and d9 are constants fixed %rlth respect to the body and u is the 

particle velocity relative to the body. With 

« 

u - T, + Tj9f (x) + n • 0 (36) 

and dr^ ■ dx/y, dtj “ fdd, Eq. (35) reduces to 

it * 7^ ‘’'‘"j ^thtri* [♦*<*>‘>^“] • ® 

where the velocity components x/y(x) and ^f(x) are functions of x, 9, and t 
as given by Eq. (29) (in which is freely replaced by 9) with appear- 
ing as a parameter. Equation (37) is thus a first-order linear partial differ- 
ential equation for pdn. 

The boundary condition required may be obtained by applying the lav of 
conservation of mass across the shock at the outer edge. Thus, the mass-flow 
rate of the free stream passing through dr^dij at (x^, 9^) is 

dij - p.(t - v^) • (-n dTjdTj) - f^If^ + (0 + y)cos 

+ A(xj - h + f^fpeoe ♦iJdxi d#^ (38) 
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( 19 ) 


wheraas tha aaaa-flow rata iMadlataly bahind the shock la 



p(Xi,^p|u|f^ dn ddi 



d#, 


Consarvatlon of aass daaaiuls ddi^ * ^^2 and tbia yields tha raquirad boundary 
condition that at x ■ x^» 4 ■ 4i* t ■ t^: 

^ pdn - + (0 + y)cos + e(x^ - h + f^fpcos 0i]dx£ (40) 

Tha solution to Eq. (37) with bouiuiary condition Eq. (40) is sought in 
tha form 

pdn - (A(x) V [0B(x) ^(x) + tD(x) + yE(x)]cos 0}dx^ (4l) 

Substituting Eq. (41) into Eqs. (37) through (40) » and equating like terms in 
0» 0» Yt Y> etc., wa obtain the following set of ordinary differential aqua- 


tions for tha dataiminatlon of A to E: 

(AA)*-0; (AA)^ - f^fj (42a) 

A 

(BA + AB)’ + ^ - 0 ; (BA + AB)^ - f^ (42b) 

A A 

(CA + AC) • + ~ + ® 0 ; (CA 4 AC)^ - f^(x^ - h + fifjl (42c) 

(DA + AD)' + “ - 0 ; (DA + AD)^ - (42d) 

A A jyf n A A 

(EA + AE)'+^ + ^- 0; (EA + AE)i-0 (42e) 


Equations (42) can be solved successively, yielding 
f(x.)f'(x ) 

A m ■■ ■■ ■ ■ 

w(Xj) 

f(x.) 

B(x,x.) ■ D(x,x.) - ^r- — - {1 - u(x.>t(x.)f'(x.)S 2 (x,x.)J 
^ J- W (x^) 111 1 


(43a) 


(43b) 
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. f(x.) (x. - h (x. - h) 

II + f"(x )m (X )1 

i I — — j f(,)d, ^ ^ 

"'J J y( 0 f *(0 J 


2£’(xJ 


P*(Xj) 


Sq(x,x^) 


"i 1 

£(x.)f*(x.) I 

^ T, — So<x»x )Sa(x.x )> 

M (Xj^) ^ ^ I 


f(x^) 




f'(x,) H + f’’(x )m'(x )) 

-5— i- (f(x) - f(x,)l j-i- i— S„(x,x.) 

1 U (Xj^) * 


w (x^) 


f(Xj)f'(Xj) f'(Xj) (•* (•* 

— j^So<x,x^)S,(x.x,)-^-— J 


mi. 

w(») 


*l -*1 

We ere now in e poeltlon to celculete the centrifugal pressure correction 
The noraal coaponent of accelerstion isi froa Eqs. (22c), (29)* and (32)* 

a^(x,^*t;Xi) - <(x)A* + 2|etc(x)AB + &(k(x)AC + A] 

■f yk(x)AD y|k(x)AE + ^ 


where 

x(x) - y*(x)f'(x) 

is tne curvature of the body surface along direction (tl^ curvature of the 
body surface along direction course* equal to -I/f(x) for the 

body of revolution). Putting Eqs. (41) and (44) into Eq. (34)* using the 
solution Eq. (43), and expanding and neglecting terns of 0(8^. dy), we 
get the centrifugal pressure correction: 


(43c) 


(43d) 


(44) 


(45) 
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X X 

/ {w(Of(OU - f'^(01 

- f*(C)f’(OSj(x.O)d5 + & cos J*[y*(OC(x.5) 

+ 2fa)f(OC(x,Ol<iC+0 cos ^ f(x) + Y cos ♦ j [n"(5)E(x,C) 

+ 2£(Of’ a)E(x,e)ldt + ; cos ♦ ^ f d? 

*o 

Adding Eq. (46) to Eq. (17), we get for the surface pressure 
?(*.♦. t) - Pq(x) + cos ♦{[0(t) + Y(t))Pi(x) 

+ 0(t)[p2(x) - hPj(x)] + Y(t)Pj(x)} 


tdiere 


- a"(6)f’(OS,(x,0 - 2f(Of 


' . r f (Odcjde 

y(Ti)f^(n) 


“ £(!)■ / ^ 'f"W / tw'(OE(x.^) + 2f(0f’(0E(x,01d5 


(46) 


(47) 


(48a) 


P^(x) - y^(x)f^(x) +f^ J[ v(C)f(Of’(OdC 

P^(x) -2y^x)f’(x) + J*{y(Of(C(l “ f’^(5)l“ f^(Of’ (C)S 2 (x.O )dC (48b) 
? 2 (x) - 2y^(x)f’(x)tx + f(x)f'(x)] +f(x)+f^ J |5f^Oy(OH * f'(0] 

+ 2y(0f"(0f’(0 - f(OSo(x,5) -h So(x,OS 2 (x,U 


(48c) 

(48d) 


The expressions in Eqs. (48) may be further simplified by elementary operations 


on the Integral terms, yielding 


Pi(x) - 2y"(x)f(x) +1 f y(0f(0[i - 2r’iO]d^ 

P 2 (x) = y^(x)[f(x) + 2xf ’ (x) + 3f(x)f’*(x)l 

**' ?(x) " Jz(x) - Js(x) + J^(x)] 


(49a) 


(49b) 


12 


(49c) 



**®^’*^ " fix} f(x} ■*“ ~ J Js(C)(f(x) - fCOJdf, 

-/> Js(e)S,(x,Odc| 

vh«re 


Ji(x) 

Jj(x) 
J,(x) 
Ji,(x) 
J j(x) 



(50) 


In aunaiary, with the shape function f(x) of a body of revolution given, 
the Newton-Busemann surface pressure in unsteady flow can be obtained by quadra 
ture. The function P^(x) in Eq. (49a), which represents the perturbation 
pressure at snail incidence in steady flow, is identical with that given by 
Hayes and Probsteln* (whtn a misprint in Ref. 2 is corrected: t’.c upper limit 

of the Integral (3.8.17) of Ref. 2 should be Xj Instead of 1). The functions 
P 2 (x) and Ps(x) are new, giving the out-of-phase pressure components arising 
from unsteady motions in 0 and y, respectively. 
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V. Stability of Oscillating Bodies of Revolution 
The general results tor unsteady surface pressure obtained in Sec. IV 
will now be used to study the dynaalc stability of bodies of revolution. An 
in Ref. 1, the principal motion variables are taken to be the angle of pitch 
6(t) and the angle of attack a(t), where a(t) - 6(t) + v(t). Equation (47) 
is then recast as 

p(x,^,t) - Pq(x) + cos ♦{o(t)Pi(x) + 6(t)[Pj(x) 

- P,(x) - hPi(x)] + i(t)P,(x)} (51) 


The normal force and pitching-moment coefficients are again defined as usual 


N M 

S " (l/2)p„U«*S^ * ^ " (1/2)p*U«ZS^£ 


(52) 


where is the base area of the body. N is the normal force, and M is the 

moment of surface-pressure force about the pivot axis C. The various stability 
derivatives are defined as in Ref. 1. Thus, 





^ p(x,i^,t)co8 dS , (|» - 0,a,a 

p(x,(fr,t) [x- h+ f (x)f’ (x) jeos <f dS , 


(53) 

e 

■ O.ci.d 


where S is the lateral surface of the body of revolution. The stability 
derivatives introduced above are also related through the axis-transfer rules 
given by Eq. (59) of Ref. 1. In this paper we shall study in detail the 
stability of the pitching motion of classes of bodies of revolution in recti- 
linear flight [y(t) ■ 0. 6(t) ■ a(t)]. In this case the stiffness derivative 

-C^p and the damplng-ln-pitch derivative take the form 

1 

-C„Q - J Pi(x)f(x)[x - h + f(x)f’(x)]dx (54a) 

"S " f’flT / ‘ hPi(x)]f(x)[x - h + f(x)f'(x)]dx (54b) 
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Circular Coneu 


% 


% 


For a cone 


f(x) - X tan T (55) 

where the conetent t Is the seadvertex angle of the cone. Froa Bq. (49), 
we get 


Pq - stn* T 
whence frMB Eq. (54) 


• 2 sin T cos T , ?2 ■ 3x tan r 


-Cmg 

-Cm| 



(56) 


(57a) 

(57b) 


The stiffness derivative formula (Eq. (57a)) according to Newtonian theory 
has been extensively demonstrated (cf. Figs. 8-11 of Ref. 3) to yield very good 
agreesent with the results of experiments for a 10* cone In air flows at Mach 
nuabers greater than 5. Equation (57b) for the damplng-ln-pltch derivative of 
a cone Is, on the other hand, new and heretofore unverified. It contains the 
centrifugal force correction due to the curved paths followed by fluid particles 
in unsteady motion and generalizes the result of Mahood and Rul** to Include the 
dependence on pivot position h. 

Although, for a cone, a centrifugal force correction Is not required for 
steady flow, such a correction i^ required for unsteady flow and, in fact, is 
fully as important as the impact pressure contribution. This is clearly demon- 
strated In Fig. 2, in which the full damping-ln-pitch coefficient with centri- 
fugal force correction Included (Eq. (57b)) Is compared with the coefficient 
that results from consideration of the Impact pressure alone, namely (from 
Eq. (17)), 
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T I 2 3 


h co«' T+ (h co®^‘ t)‘ 


The difference between the two results In Fig. 2 Is e direct eessure of the 
centrifugal force contribution, which, analytically, takes the form 




centrif 


_ri.i 

T L4 3 


h cos* T 


It will be seen in Fig. 2 that the centrlfugal>-force contribution (Eq. (59)) 
can be of the sasMfi order of magnitude as the impact pressure contribution 
(Eq. (58)). This result casts doubt on the validity of various unsteady 
Newtonian flow theories that are based on the Newtonian impact pressure alone. 
For example, Ericsson's theory,^ in which -Cg, is calculated from the 

q 

Newtonian liq>act pressure alone and then added to an estimate of to 

& 

give the equivalent -0^ , must be Incorrect In at least the separate coa- 

8 

ponents, if not the sum, since both components contain centrifugal-force 

contributions. These separate contributions can be deduced from the results 

presented in Eq. (49). Thus, the centrifugal-force component in -C,^ (which. 

a 

for this coefficient, is also the total since there is no impact pressure com- 
ponent) is calculated from P,(x) alone, and takes the form 




centrif 




which is seen to be exactly one half of the total centrifugal-force component 

in -Cm. (Eq. (59)). It may be easily verified that the centrifugal-force 

6 

component in -C_ contributes the other half. This result, which is also 

q 

true for the wedge airfoil (cf. Eq. (51) of Ref. 1) may be summarized as 
follows: the centrifugal force components in -C^ and -C^ are equal (for 

q a 

cones and wedge airfoils) , and therefore each coefficient contributes exactly 
one half of the centrifugal-force coeqtonent in -C^. . 
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There It yet another aspect of this result that warrants consideration. 

The fact that has a significant centrifugal-force coRqmnent means that 

even in the Newtonian lladt, an accurate estimate of -CL cannot be derived 

A 

on the basis of a local analysis, such at laquict theory, or In fact any theory 
In which local pressure Is said to depend solely on the local flow velocity. 

On the contrary, the presence of a centrifugal-force term is an indication of 
the Influence of flow history on the local pressure. Since a local analysis 
for -C^ is Inadequate even in the Newtonian limit, this should give pause 

q 

to those proposing to predict aerodynamic stability derivatives at lower speeds 

on the basis of clmllar local analyses. 

Also plotted in Fig. 2 are the results of measurasents for the damplng- 

In-pltch coefficient carried out In ulr flow at a free-stream Mach number of 

10 by Hobbs^ and the results of ntaaerical computations based on an Invlscld 

gas-dynamic theory carried out by Brong.* It will be seen that results based 

on the complete Newton-Busemann theory are in excellent agreement with the 

experimental results for the highest Reynolds number presented; the numerical 

results of Brong agree better with the experimental results for lower Reynolds 

numbers. We believe that it Is partly fortuitous, and not an Indication of 

the superiority of Newton-Busemann theory, that results from the latter theory 

should be In closer agreement with the experimental results at the highest 

Reynolds number than those of the gas-dynamic theory. In Ref. 8, dealing with 

gas-dynamic theory. It has been shown that both Increasing the free-stream 

Mach number Mm and decreasing the ratio of specific heats y tend to increase 

the damplng-ln-pitch coefficient -CL.. Consequently, Newton-Busemann theory, 

6 

the theory that results from taking the double limit Mm Y 1 of gas- 

dynamic theory, should overestimate -C_. for finite values of M„ and for 

“e ^ 

values of y greater than unity (y ■ 1.4 for air). Physically, this is 
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because In the Newtonian Italt the ehock Inycr becoace inflnltcelmalty thin* 

thue elimlnatlx^ one aechanlaa of dynwRlc dcatabtllsatjon, nanely, tho pli.tHo 

lag reaultlng froa preaaure algnala traveling to and reflecting froa the bow 

ehock. Since NefwtonrBueeaann theory ehould have overpredicted -Cg|.» the 

6 

close agreement evidenced in Fig. 2 between results froa Ne%rton-Buaaaann theory 

and the experlaental results at the highest Reynolds number must be attributed 

to the presence of the boundary layer on the body In tbs experiments, which, 

in the usual way, may be considered to have formed an effective displacaaent 

surface. The effective thickness of the cone would have increased thereby, 

providing a mechanlsa for a correspoxuling increase in **C_ . Thus, the 

6 

Newton>Busemann theory may have overpredicted -C^. by Just the right amount 
to account for viscous effects. 

Finally, we note in Eq. (57b) that analogous to the wedge airfoil (Eq. (62)), 
Ref. 1), the damping-in-pltch coefficient -Cgi^ for a cone is always positive 
according to Newton-BusMsann theory. Hence, the pitching motion of a cone of 
any thickness and about any axis position is dynamically stable in the 
Newtonian limit. 

Bltmt Cones 

To investigate the effects of nose bluntness on the dynamic stability, we 
consider a class of bodies in which, for each body, the nose consists of a 
npherlcal cap, and the afterbody is the frustum of a circular cone which joins 
the spherical cap smoothly and with a comson tangent (Fig. 3) . Let the radius 
of the spherical cap be Rfj and the radius of the base of the cone be R^. 

Then the class of bodies under consideration is described by the following 
parameters: the saslangle of the ue t, the ratio A ■ 

(actual) length of the body (from nose to base) 1 ■ 1. Let 
s • (1 - sin T)/ain t; then 
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X tan T 
1 > flX can T 


» 


Rn • 


Lat J “ aR)|; then the length c of the full cone la c • 1 + J and 
d)tan t. The eqtiatlon of the body la given by 

1 (2R|^x - X*)'/* ; 0 < X < j sin t 

(x j)tan T ; j sin t < x 5 1 


( 61 ) 


Each of the etability derivatives consists of two parts, the first deriving 
from the pressure on the spherical surface ( ),. and the second from the 
pressure on the conical surface ( )^. Thus. 

(■®"0 ■ (■'»«). 


(62) 


On Che spherical surface, using Eq. (61a). t#e get from Eq. (49a). 

Pj(x) - I (I - xj)(2xi - ^0 5 xi s ^ 5 I - sin (63) 


It can also be shown, either by physical reasoning based on Eq. (47) or by 
direct coBq>utatlon in Eq. (48), that 

P2(x) - R„Pi(x) (64) 


Reasoning physically, we note in Eq. (47) that the contribution due to the 

s 

nose cap to the pressure component proportional to 6 must be zero when the 

axis h. about which the body turns, passes through the center of the spherical 

nose cap (i.e.. h “ for then the curve defining the nose-cap surface turns 

on Itself and there la no change in the body geometry, so far as the invlscld 

external flow is concerned, to alert it that the nose cap Is turning. We sec 

from Eq. (47) that the pressure component proportional to 6 will be zero at 

h - Rg if and only if Eq. (64) is true. Using Eqs. (63) and (64) to evaluate 

the contribution of the nose cap to -G, and -C^. , we get 

6 6 
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■ 3 X* CO*** t(Rjj - h) I 


(65) 


(-C«§)^ -f X* cos'» t(R^ - h)*j 

On the conical afterbody, k(x) " 0, hence P^(x) - 2 eln t coa t and 

PjCx) - (3x + j(l + 2 sin* T)]tan t. Theae ahow that, owing to the vaniahing 

of K» the pressures on the conical afterbody depend on the local position x, 

but not on the body shape upstream of x.^ Making use of this property in 

calculating the contributions to •Cn and -Cn. due to the conical afterbody, 

0 0 

we are permitted to build up the contributions of the afterbody out of the 
contributions fr<» full cones; that is, those of a large cone with origin at 
the apex of the conical frustum and base at the base of the frustum (cf. Fig. 3), 
less those of a small cone with origin at the same apex and base at the foot 
of the frustum. We get 


- a + J)[n»o) - 


( 66 ) 


where 6 « X cos t and Xq ■ (h + j)/(l + j) is the pivot position from the 
apex of the large cone relative to its length c (Fig. 3). The functions 
F and G are defined in Eqs. (57). It should be apparent that each of the 
first terms in Eqs. (66) represents the contribution of the large cone, and 
that each of the second terms represents the (subtractive) contribution of the 
small cone. 


^Thls property, which holds for conical bodies only in the Newtonian limit, 
makes the theory a local one in the sense that the pressure at a point x is 
"unaware" that the frustum is not a full cone. 
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In the greeter pert of the literature dealing with the atablllty of blunt 

conee»^*^ the pltching-noaent coefficient end the atability dertvatlvca are 

defined with the length c of the extended sharp cone being used as the 

cherecteristic length rather then the actual body length that is used in 

this paper. We note here that the stiffness end deaping- in-pitch coefficients 

fron Refs. 7 end 9» formed with c as cherecteristic length, are equal to 

-(1 + j)”^Cn , -(I + , respectively, in our notation. The effects of 

9 6 

bluntness on the stability derivatives may be exaained by forming the ratios 
of the derivatives with their counterparts for sharp cones. Combining Eqs. (65) 
and (66), we get 


fa + ^ hw'' 3 “ *0 T)J 

I Jr^.0 


> (67) 


<1 + 

fa + 


- 1 + 


G(x. 


? 5 (6 - *o CO." T)*j 


*N-« 


For small bluntness, X • R{)/Rg « 1, hence 6 « 1, and the right-hand sides of 
Eqs. (67) both reduce to 1 *1- O(X^), shoving that the effect of nose bluntness 
should be negligibly mall for small bluntness. This theoretical prediction is 
borne out by a coaqiarison of experimental results with coBq>utations based on 
Nevton-Busemann theory for blunted cones (with t fixed at iO* and x^ fixed 
at 0.6). The co^>arisons are presented in Figs. 4 and 5, where the experimental 
points have been taken from an extensive compilation provided by Khalid and 

East in Ref. 9. We note in Fig. 4 that the experimental results for C_ 

®() 

would appear to confirm the initial behavior predicted by Newton-Busemann theory 
for values of X up to about 0.07. Similarly, in Fig. 5, as predicted by 
Ne%rton-Busraann theory, experlatental results for initially do not vary 
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significantly with X for values of the paratMter up to about 0.05. As the 
bluntness paroteter Increases, however, Newton-Busesuinn theory is Incapable of 
predicting the abrupt and large variation with bluntness that is observed 
experlnentally in air flow with y ■ 1.4 for both the stiffness and the 
daiBping~ln>pltch coefficients. This incapacity is synptoaatlc of Newtonian 
flow theory, both steady and unsteady, in cases where a principal requirenent 
of the theory is violated, nasMly, that the shock layer renain everyidtere thin. 

As seen in Fig. 6 for spherically blunted cones, as the nose bluntness increases, 
the shock over< xpands around the nose cap, creating a shock layer of significant 
thickness on the conical afterbody. Possible amendnents to the theory are dis- 
cussed in Sec. VI. It is of Interest to note that coiq>utatlons based on a com- 
plete gas-dynanlc theory^'' — results from which are also plotted on Figs. 4 
and 5 — suggest that gas-dynamic theory is capable of capturing the major fea- 
tures of the dependence on bluntness. Curiously, henirever, as computed, the 
numerical results miss the initial behavior of the derivatives as a function of 
bluntness parameter. 

Slender Powr-Law Bodice 

For slender bodies, f(x) ■ 0(b) with b « 1. When terms of O(b’) are 
neglected compared with terms retained of 0(b) , the formulas for surface 
pressure (cf. Eqs. (49) and (50)) simplify to 

Po(x) - f’*(x) + J f(x)f'(x) 

Pj(x) - 2f (X) + I f(Odt 

° , ( 68 ) 

Pj(x) - f(x) + 2xf’(x) + J [2KHO 

- 2 f f(n)dn - f f*(n)dn)d5 

*o *o 
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For po%fcr-’law bodioa of revolution 

f(x) - bx" 


(69) 




* 


we ^et 

p^(,) . b* ailHjzJI 

?3l(x) - 

Pj(x) “ bMjX*^ 

where 

., » non * 1) 

"» 3 n + 1 

Mj - (2n + 1)^1 3(a + l)tn^-f 2)(3n\ li 

The d«Bplng-in>pitch derivative ia 


(70) 


(71) 


-<=•6 ■ ir^ - * *llr^ ^ ^ «”> 

For Newton-Buaeaann flow tlMory to apply, it ia required that PqCx) a 0, hence 
n > 1/3. The damping-in-pitch derivative Eq. (72) haa a ainiatuB value given by 

which is plotted in Fig. 7 aa a function of the power n. It ia acen t'lat 

increaaing the power n, that ia, making the ahape of the body aore concave, 

tenda to decrease (~C..) . In thia sense it aay be aaid that increaaing the 

6 Bin 

convexity of alender power-law bodiea of revolution haa a favorable effect on 
the dynamic stability of the bodies. This conclusion is aiBilar to that 
reached in Ref. 1 for el«;der power- law airfoils. However, for bodies of revolu- 
tion, (-CL*) > 0 for all n. Implying that the pitching motions of alender 

6 Bin 

power-law bodiea of revolution are always dyn«ically stable. In contrast. 
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slondor power-law nlrrullN of concava aha pc nay be dynamically unalable if 
the nurfarc curvature cxceedn o certain limit. 

VI . Concluding Rmmarka 

In thia pmpnr, wm hava davalopad a complata unacaady Mawton-Buaamann flow 
theory for oaclllatlng bodiaa of revolution of arbitrary ahape having aharp or 
blm:^ aoaea. Exact forwilaa in cloaed fora are given for the atability deriva- 
tivea of thcae bodlea; they require only nuawricel quedraturea of terms involv- 
ing the body shape. 

The centrifugal force arising from the curved trajectories followed by the 
fltiid particles in unsteady flow is alunm to be very important. Its contribution 
cazmot be neglected, mvKtt for a circular cone for which the centrifugal force 
contribution la sero In steady flow. With this correction included, the theory 
ia shown to be in excellent agreement with experimental results for sharp cones 
and for cones with small nose bluntneas. 

However, as expected from the Incapability of steady Newton-Buaemann flow 
theory to predict the surface pressure on blunt-nosed bodies in air flow (with 
Y - 1.4), the unsteady Newton-Buaemann flow theory also gives poor results for 
stability derivatives of bodies with moderate or large nose bluntncss. Thla Is 
primarily due to the violation of ' the principal requlreawnt of the theory in 
Much flow, namely, that the shock layer remain everywhere infinitesimally thin. 
Improvements nay result by us''Tig the thin-shock- layer theory. In tihlch solutions 
to the full gas-dynamic equations are sought in the form of a power series in 
the thickness of the shock layer. The unsteady Newton-Buaemann flow given here 
would then provide ttM leading term In such a rational expansion scheme. For 
configurations such that Che shock layer remains everywhere thin, the chin- 
shock-layer theory may be expected to give only mall corrections to the Newton- 
Busemann theory. Such Is Che situation In the case of an oscillating delta wing. 
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^ Ind««d» It can b« ahoim^^ chat c)m thln-ahock-> layer theory Tor oHcillatlnK <lolt.i 
vlngs*^ givaa raaulta nlaoat identical with those obtained with unateady Newton- 
Buaaaumn theory. 

An alternative but eaplrlcal approach* which any produce inproveaent in 
predictive potser over the preamit theory* would be to apply the cmicept of 
eabedded Newtonian flow,^*^* The application ia aade feaaib.le by the avail- 
ability now of a coaplete deacription of the r nateady flow field baaed on the 
Newton-Buaeaaim concept. 
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FIGURE CAPTIONS 


Fig. 1 Oscillating body of revolution. 

Fig. 2 Comparison of theory and experiment for daaplng-ln-pitch derivative of 
a 10* sharp cone. 

Fig. 3 Geometry of blunt cone. 

Fig. 4 Variation of 10* cone stiffness derivative with bluntness — comparison 
of theory and experiment. 

Fig. 5 Variation of dan^ing-in-pitcb derivative with nose bluntness — 
comparison of theory and experiment. 

Fig. 6 Variation of shock- layer thickness with bluntness of a 12.5* cone in 
hypersonic air flow. (Courtesy of G. Malcolm, NASA, Ames Research Center.) 

Fig. 7 Variation with power n of minimum damplng-ln-pltch derivative for 
slender power-law body. 
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